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ABSTRACT 

Recent surveys have revealed a lack of close-in planets around evolved stars 
more massive than 1.2 M & . Such planets are common around solar-mass stars. 
We have calculated the orbital evolution of planets around stars with a range of 
initial masses, and have shown how planetary orbits are affected by the evolution 
of the stars all the way to the tip of the Red Giant Branch (RGB). We find that 
tidal interaction can lead to the engulfment of close-in planets by evolved stars. 
The engulfment is more efficient for more-massive planets and less-massive stars. 
These results may explain the observed semi-major axis distribution of planets 
around evolved stars with masses larger than 1.5 M & . 

Our results also suggest that massive planets may form more efficiently around 
intermediate-mass stars. 

Subject headings: Planetary systems — Stars: evolution 



1. INTRODUCTION 



Observationally, due mostly to the inapplicability of high-precision Doppler techniques 
to stars with spectral types earlier than late-F (which have large rotational velocities and a 
small number of spectral lines), little was known until recently about the frequency of planets 
around the more-massive stars. This situation has now changed, as survey s have been ex- 



tended to searches for planets around more-massive stars in an evolved stage (IDollinger et al. 
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One of the most important trends that these su rveys have revealed is the lack of close 



in planets orbiting stars with masses M > 1.5 M (j Johnson et al.ll2007bl : ISato et al.l 12008 
Wright et al.ll2009l ) despite the fact that these planets are found around ~20 % of the Main 



Sequence (MS) stars with M < 1.2 M^. The frequency of planets seems al so to be higher 
around intermediate- mass stars (ILovis fc Mayorl 120071 : 1 Johnson et al.ll2007bl ). Furthermore, 
contrary to the correlation between metallicity and probabil i ty of planet-hostin g found for 



solar- mass, main-sequence stars (e.g.. iFischer &: Valentill2005l : ISantos et al. 



massive planet-hosting stars do not exhibit higher metallicities (e.g.. iPasquini et al.l 120071 ) 



200 5]), the mor e 



On the theoretical side, it has been shown that the form a tion of Jupiter-m ass planets 
around M stars may be hindered (e.g., lLaughlin et al.l 12004 ; llda fc Lml 120051 ). while the 



probability that a give n star has at least one g as giant increases linearly with the stellar 
mass up of 3 M & (e.g., Kennedy fc Kenyonl 120081 ) . 



One possibility is that the observed difference in the orbital distribution of planets found 
around intermediate and solar-mass primaries is due to the evolution of the star. Since all 
of the planets orbiting red giants or subgiants with M > 1.5 M have semi-major axes 
a > 0.5 AU it has been suggested that the pla nets might be engulfed as the star evolves off 



the MS ( I Johnson et al.ll2007bl ; ISato et al.ll2008l ). This possibility is, however, often dismissed 



(Johnson et al. 


2007b. 


2008; 


Currie 


2009) 



the observed differences in orbital distribution are primordial, and and they are a consequence 
of the pla net formation mechanism around more massive stars. Along these lines, it has been 



shown by lCurriei (120091 ) that the dependence of the lifetime of the gaseous disk on the stellar 
mass could result in halting the inward migration of planets around high-mass stars, thus 
explaining the observed lack of short period planets around these stars. 

The point we are making in the present paper is that before stellar evolution can be ruled 
out as the mechanism behind the observed semi-major axis distribution of planets around 
evolved stars, detailed modeling of the orbital evolution needs to be performed. In other 
words, in order to determine the potential role of the stellar mass in the planet formation 
process, the effects of the evolution of the star on the observed orbit distribution around 
giants have to be correctly isolated. This is precisely the goal of this paper. 
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2. THE EQUATIONS 



There are several competing processes that affect the orbital distance between the star 
and the planet as the star evolves off the MS: the changes in the mass of both the planet 
and the star, the gravitational and frictional drag, and the tidal force. 

To determine the rate of change in the planet's mass, we consider a planet of mass M p 
and radius R p moving with a velocity, v, in a circular orbit (e = 0) around a star of mass 
M*. Since the planet is moving supersonically through the matter ejected by the giant star, 
it accretes mass. The accreti on rate onto the planet, M p ^ is given approxim ately by the 



Bondi-Hoyle expression (e.g.. lBondi fc Hoyldll944l ; iRuderman fc Spiegellll97ll ). 



M, 



p | acc 



nR A pv 



(1) 



where p is the density of the environment and Ra is the accretion radius (Ra = 2GM p /v 2 
with G the gravitational constant). At very short distances (where Ra < R p ) we have 
replaced R\ by RaR p to correct the geometrical radius by gravitational focusing effects. 

At the same time, the planet's surface is being heated by radiation arising from the 
shock front and from the stellar surface. This heating can lead to evaporation of surface 
material. We estimate the evaporation rate, M p i ev as in Villaver & Livio (2007; Eq. 9). 
The temperature at the planet's surface has been taken to be the maximum between the 
radiative equilibrium temperature of the planet (see e.g., Eq. 5 in IVillaver fc Livid 120071 ) 
and the temperature of the shocked gas T s h = (3m#/16A;{,)t> 2 (estimated from the Rankine- 
Hugoniot conditions for an adiabatic shock where is Boltzman's constant and the mass 
of the hydrogen atom). 



The rate of change in the planet mass is thus given by 



{M p , acc - M p 

CV J 



(2) 



The rate of change of the stellar mass is simpl y M* = — M m ] r , where M m i r is the stellar 
mass-loss rate. Using Reimers' law for Red Giants (lReimerslll975l ). 

L^R 



M, 



mil 



4x W~ 13 r) R 



[M Q yr- 1 ] 



(3) 



where L*, i?* and M* are the stellar luminosity, radius, and mass respectively (in solar units) 
and t)r is the Reimers parameter, which we take throughout this work to be t^r = 0.6 

Conservation of angular momen tum gives the equation for the rate of cha nge in the 
orbital radius of the planet (see, e.g., Alexander et al.lll976t iLivio fc Sokerlll984l ) . 



M* + M p 
M, ■ M,, 



[Ft + Fg 



(4) 



-4- 



where Ff and _F 5 are respectively the frictional and gravitational drag forces and (a/a)t is 
the rate of orbital decay due to the tidal interaction. 

The gravitational drag force, F g , arises from the eddying motions that are set up in 
the fluid by the passage of the planet. It is a consequence of the gravi tational intera ction 
of the planet with a gaseous medium. The drag force is given by (e.g., lOstrikerl Il999l . and 
references therein) 



47T- 



- P I 



(5) 



where / is a time-dependent function of the Mach number. The numerical results of lOstriker 
(119991 ) show that for the Mach numbers encountered here, / is approximately constant and 
has the value I ~ 0.5. 

The loss of angular momentum associated to the frictional force Ff is proportional to 
the surface area of the planet exposed to the flow and it can be expressed in the form (e.g., 
Rosenheadl Il963l ) 

F f = ~C dP v\R 2 p , (6) 

where C d ~ 0.9 is the dimensionless drag coefficient for a sphere. 

Finally, the angular momentum loss associated to the tidal term (a/a) t arises from the 
additional force (besides the gravitational pull between the two centers of mass) resulting 
from the non-spherical part of the mass distribution from the tidally distorted companion. 



Zahn 


1966. 


1977. 


1989) 



produce tidal friction is turbulent viscosity (e.g., 
timescale is determined by the effective eddy viscosity, with eddy velocities and length scales 
given appro ximately by standard mixing length theory if convectio n transports most of the 



energy flux (jZahnlll989l ; IVerbunt fc Phinneylll995l ; iRasio et al.lll996l ). The tidal term is given 
by 



Td M* V a 



(7) 



with M env being the mass in the convective envelope, q = M p /M*. and Td the eddy turnover 
timescale, given in the case of a convective envelope (iRasio et al.lll996l ). 



Td 



M env (R* — R 
3L* 



1/3 



where R env is the radius at the base of the convective envelope. The term / in Eq. ([71) 
is a numerical factor obtai ned from inte grating the viscous dissipation of the tidal energy 
across the convective zone. IZahnl (1l989l ) used / = 1.01(a/2) where a is the mixing length 
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parameter. IVerbunt fc Phinneyl (119951 ) confirmed that observations are consistent with / «1 
(specifically, they obtained values in the range 0.5 < / < 2 with a preference to values in the 
upper half of the range) as long as r,; < P with P being the orbital period. We therefore 
used / = (P/2r^) 2 to account only for the convective cells that can contribute to viscosity 
when Td > P/2, otherwise we take / = 1. 

The evolution of orbital eccentricities is beyond the scope of the present paper. The 
damping of eccentricities due to t idal forces may eventually lead to a nar row distribution of 
e- values for close-in planets (e.g., iRasio et al.lll996l ; I Jackson et al.ll2008al ). Note that the 



initial value of the eccentricity has little effect on the orbital decay rate (e.g., 
2008bh . 



Jackson et al. 



Calculation of the tidal term requires a knowledge of the structure of the star. We have 
used detailed stellar models provided to us by Lionel S iess. These were calculated based 
on the stellar evolution code STAREVOL described in ISiesd (120061 ). We have used stars 
with MS masses of 1, 2, 3, and 5 M , solar metallicity, and a mass-loss prescription with a 
Reimers parameter of T}r = 0.6. 



3. RESULTS 

To determine the evolution of the planet's orbit we integrate Eq. H] along the RGB. The 
stellar evolution timescales (in the absence of significant mass loss) are set by the rate of 
consumption of the nuclear fuel. Since the nuclear burning has been exhausted in the core, 
during the RGB hydrogen burning continues in a shell outside the helium core, which now, 
devoid of energy sources, is contracting and heating up. As the core contracts the envelope 
expands and cools. The stellar effective temperature decreases while the star's radius and 
luminosity increase. We have updated the stellar parameters at each time step of the orbital 
evolution. 

As the star evolves off the MS the wind velocity decreases (Kdnd ~ 5-10 kms -1 ). The 
termination shock which marks the position where the stellar wind interacts with the ISM 
cools down and the temperature of the medium at the planet's locat ion has an upper limit 



provided by the stellar effective temperature (e.g. JVillaver et al.ll2002l ). We used T = 2500 K 
as the temperature of the ambient medium. The Mach number M = v/c s is then always 
M > 3 for the orbital distances relevant to this work. 



The density of the environment has been computed as 

M m i r 
P 4vra 2 K v ind 



(9) 
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unless M m i r = 0, in which case we took a number density of p = 0.001 cm 3 . We use a wind 
velocity of 5 kms _1 along the RGB. 

We have computed the orbital evolution for a range of initial orbital distances and planet 
and stellar masses in order to determine the minimum initial orbital distance for which a 
planet will avoid being tidally captured by the expanding star. For illustration, some of the 
computed orbits (dash-dotted lines) are plotted in Fig. CD, where we have used a planet with 
a mass of M p = 1 Mj orbiting stars with different MS masses to the end of the RGB phase. 
The evolution of the stellar radius is shown as a solid line. The panels, from top to bottom 
and left to right, are for stars with MS masses of 1, 2, 3, and 5 M & respectively. Note that 
the scale on the axes is different for each panel. 

At large distances from the star, the densities involved are low, and the drag terms 
associated with the forces Ff and F 9 in Eq. (J3J) play a negligible role in the evolution of the 
orbit. Moreover, since the accretion rate onto the planet is always small compared to the 
stellar mass- loss rate, the first term in Eq. (jlj) is dominated by M*. The temporal behavior 
of the orbit is then mostly governed by the relative importance of the terms associated with 
the stellar mass-loss M* and the tidal interaction (a/a) t . 

Red giant mass-loss rates are somewhat uncertain. As noted in §2 we have in the present 
work estimated the mass-loss rate by the Reimers prescription with t]r = 0.6. This seems 
to reproduce fairly well the observations of individual RGB stars. The peak RGB mass-loss 
rates are higher for lower-mass stars (10~ 8 for 1 M versus 10~ 10 for the 3 M & ), and the 
lowest mass stars also reach the largest radius at the tip of the RGB. 

Figure [1] demonstrates the three possible outcomes of orbital evolution: (i) Beyond a 
certain initial orbital separation, the orbital separation simply increases, due to systemic 
mass loss, (ii) There is a range of initial orbital separations for which the orbit decays, but 
the planet avoids being engulfed, (iii) Inward from some critical, initial orbital separation, 
the planet is engulfed mostly due to tidal interaction. 

With the purpose of quantifying the influence of the planet's mass, we have also inte- 
grated the orbit of planets with masses of 3 and 5 Mj (see Figs. [2] and [3]). The top panel of 
Fig. [2] shows the RGB evolution of a 2 M star along the HR diagram (Fig. [3] is the same, 
but for a 3 M & star). The bottom panel shows the orbital evolution for planets with masses 
3 Mj (dotted line) and 5 Mj (dash-dotted line) together with the evolution of the stellar 
radius (solid line). Note that the bottom panel does not show the maximum extent of the 
stellar radius, but a zoom of a section at the base of RGB, to better appreciate the details. 
We selected small initial orbits to identify at which points during the RGB these planets are 
swallowed by their stars (marked by the location of the arrows). 
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For all the initial orbits that satisfy the condition a Q < i?" iax the planet gets engulfed 
by the star at same point before the end of the RGB phase. The more massive the planet 
the stronger is the tidal interaction with the star, and therefore the sooner the orbit decays 
to meet the stellar radius (see Figs. [2] and [3]). 

Some of our findings are summarized in Table [TJ where we list the minimum initial 
orbital distance for which a planet with a given mass avoids being engulfed by the star. The 
second column gives the maximum radius reached by the star on the RGB, i?™ ax (in AU), 
and the following columns give the minimum orbital distances (in AU) at which planets with 
masses of 1, 3, and 5 Mj (respectively) avoid being engulfed. For the stellar masses of 1 and 
2 M , we calculated a grid using initial orbits at steps of 0.1 AU, for the 3 and 5 M Q stars 
we used initial orbits at steps of 0.01 and 0.05 AU, respectively. 

Two important conclusions can be extracted from Table [TJ First, the tidal "capture" 
radius increases with the planet's mass; a Jupiter-mass planet is captured by a 2 M star if 
it starts at an initial orbit of 2.1 AU while a 5 Mj planet will be engulfed by the star if it 
has an initial orbit a Q < 2.5 AU. Second, the tidal capture radius decreases with increasing 
stellar mass; for a 5 Mj planet the initial orbit has to be larger than a Q ~ 3 x i?™ ax to avoid 
tidal capture around a 2 M & star, while it has to be larger than a Q ~ 2 x i?™ ax to avoid tidal 
capture around a 5 M & star. 



4. DISCUSSION AND CONCLUSIONS 

Our main goal has been to determine whether stellar evolution could explain the ob- 
served distribution of the semi-major axes of planetary orbits around evolved stars (i.e., 
semi-major axis >0.5 AU). We found that when the details of the orbital evolution are 
accurately calculated, tidal interactions constitute a quite powerful mechanism, capable of 
capturing close-in planets into the envelope of evolved stars. 

To date, there are ~20 exoplanets discovered around giant stars with M > 1.5 M . The 
host stars have radii in the range 0.02 < R* < 0.1 AU. Our models are consistent with the 
existence of the se planets at the point in the RGB evolution at which they are observed (see 



e.g., Table 8 in ISato et al.ll2008l and our Figs. [2] and [3]). However, we do not expect to find 
massive planets with a < 0.4 AU around a 2 M star with aR, > 0.1 AU (or 24 R & ). Our 
calculations provide the minimum orbital radius inside of which planets will be engulfed by 
the star at the end of the RGB evolution. 



We find that the evolution of the star alone can quantitatively explain the observed 
lack of close-in planets around evolved stars even allowing for the uncertainties associated 
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with mechanisms such as mass loss along the RGB or tidal-interaction theory. A mechanism 
such as the one invoked by Currie (2009; i.e., a lifetime stellar- mass dependency of the gas 
in the planet-forming disk that can halt migration) is not needed, although it might still be 
present. 

We find that given an initial distance at which tidal capture is possible, the more 
massive the planet, the earlier it will be captured by the RGB envelo pe. Observationally, it 
appears that giant st ars host more massive planets than MS stars (e.g. . I Johnson et al.ll2007bl ; 
Lovis fc Mayor! 120071 ) . Since we find that more massive planets are expected to be engulfed 
earlier in the RGB evolution, the fact that they are more frequently observed may point 
towards a planet-formation mechanism that favors the formation of more-massive planets 
around intermediate-mass stars. This would be consistent with a scenario in which the disk 
mas s scales with the stellar ma ss, and more massive disks produce more massive planets (see 
also Kennedy fc Kenyonl 120081 ) . 



Along similar lines, since we find a high probability of tidal capture of the planet 
by evolved stars, the higher frequency of pl anets observed around intermediate- mass stars 
( jLovis fc Mayor! 120071 : iJohnson et al.l l2007bl ) seems to imply that the efficiency of planet 
formation must be considerably higher for more massive stars, compared to their solar anal- 
ogous. 



Assef et al.l (120091 ) estimated the probability of detecting transits of planetary compan- 



ions to giant stars to be > 10 % for several of the known systems. Since tidal oribital decay 
decreases the initial orbital distance, it may increase the probability for the planet to be 
observed in transit. 

Our results sugg est that many planets may be accreted by their host star (see also 
Siess fc Liviol Il999aj jbl) . Although so far the results are based on a fairly limited sample, it 
appears that gia nt stars hosting plane ts have the same metallicity distribution as giant stars 
without planets (jPasquini et al.l 120071 ) . On the other hand, we should note that if a giant 



star engulfs a close-in planet early in the RGB, it will appear to be a giant star without 
a planet. This mechanism might perhaps provide for a partial explanation for the lack of 
correlation with metallicity. 



Finally, stars that have accreted planets may show a 
orbital angular momentum is tra nsferred to the sta r (e.g. 



ligher spin rate, as the planet's 



Livio & Sok er 2002). This phe 



nomenon has been investigated by IMassarottil (120081 ) and lCarlberg et al.l (120091 ). where they 
estimate the probability of finding rapid rotators among evolved stars. Our quantification 
of the tidal capture radius may help refine these calculations. 
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We are very grateful to Lionel Siess, who used his stellar evolution code to calculate the 
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Table 1. Minimum Orbital Radius to Avoid Tidal Capture 







M p = Mj 


&min [AU] 

M p = 3 Mj 


M p = 5 Mj 


1M 


1.10 


3.00 


3.40 


3.70 


2 M 


0.84 


2.10 


2.40 


2.50 


3 M 


0.14 


0.18 


0.23 


0.25 


5 M 


0.31 


0.45 


0.55 


0.60 
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Fig. 1. — Evolution of the orbital separation of a planet with M p = Mj (dash-dotted line) 
and the radius of the star (solid line) along the RGB. The different panels represent the stellar 
MS masses considered (1, 2, 3 and 5 M from left to right and top to bottom respectively) 
and are marked at the bottom left corner of each plot. The different dashed lines are the 
orbits for different initial orbital distances. 



-14- 




Fig. 2. — Top: The evolution of a 2 M & star during the RGB along the HR diagram. 
Bottom: The evolution of the orbital separation of a planet with M v = 3 Mj (dotted line) 
and M p = 5 Mj (dash-dotted line) for four different initial orbital separations. The evolution 
of the stellar radius along the RGB is also shown as a solid line. The arrows show the location 
at which the M p = 5 Mj planet enters the stellar envelope for the different initial orbits. 
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Fig. 3. — Same as Fig. [2] but for a star with a MS mass of 3 M & . In the bottom panel the 
star reaches the RGB maximum radius. 



